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ABSTRACT 

In the framework of Quark-Gluon String Model we calculate the inclusive spectra of 
secondaries produced in heavy ion collisions at intermediate (CERN SPS) and at much 
higher (RHIC) energies. We demonstrate that the mechanism of secondary production 
changed drastically in the energy interval \/s = 20—60 GeV and that is in agreement with 
qualitative estimates of Glauber-Gribov theory. The results of numerical calculations at 
intermediate energies are in reasonable agreement with the data. At RHIC energies 
numerically large inelastic screening corrections should be accounted for in calculations. 
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1 Introduction 



The Quark-Gluon String Model (QGSM) and the Dual Parton Model (DPM) are based 
on the Dual Topological Unitarization (DTU) and describe quite reasonably many fea- 
tures of high energy production processes, including the inclusive spectra of different sec- 
ondary hadrons, their multiplicities and multiplicity distributions, etc., both in hadron- 
nucleon and hadron-nucleus collisions [I]- [7]. High energy interactions are considered as 
proceeding via the exchange of one or several Pomerons and all elastic and inelastic pro- 
cesses result from cutting through or between Pomerons [8] . Inclusive spectra of hadrons 
are related to the corresponding fragmentation functions of quarks and diquarks, which 
are constructed using the Reggeon counting rules 

In the case of interaction with nuclear target the Multiple Scattering Theory (Gribov- 
Glauber Theory) is used and it allows to consider the interaction with nuclear target as 
the superposition of interactions with different numbers of target nucleons. 

In the case of heavy ion collisions the Multiple Scattering Theory also allows to 
consider this interaction as the superposition of separate nucleon-nucleon interactions. 
However, in this case there is no possibility to sum up all the diagrams in a rather simple 
form. The first simple classes of diagrams can be accounted as the simple expressions 
[TOl [TT] . The situation with more complicate diagrams is not so clear [T2] |14j. 

In this paper we present the QGSM results for the calculation of the inclusive spec- 
tra of secondaries produced in heavy ion collisions both at intermediate (CERN SPS, 
^snn = 17 GeV) and much higher (RHIC, ^snn = 60 — 200 GeV) energies. The 
data at GSI and CERN SPS energies [H]-[T8] are in reasonable agreement with Multiple 
Scattering Theory and with QGSM (see more detailed discussion in Section 3). However 
the similar calculations overestimate the data at RHIC energies by about a factor of two, 
as shown in Section 4. This effect was explained in [Tj5] as large contribution of inelastic 
screening correction namely in the case of heavy ion collisions. We estimate the energy 
region where these correction increase from several percents to twice as ^snn = 20-60 
GeV. 

2 Inclusive spectra of secondary hadrons in the 
Quark-Gluon String Model 

For the quantitative predictions we need a model for multiparticle production and we 
will use the QGSM for the numerical calculations presented below. 

As mentioned above, high energy hadron-nucleon and hadron-nucleus interactions 
are considered in the QGSM as proceeding via the exchange of one or several Pomerons. 
Each Pomeron corresponds to a cylinder diagram, see Fig. la, and thus, when cutting a 
Pomeron, two showers of secondaries are produced, as it is shown in Fig. lb. The inclusive 
spectrum of secondaries is determined by the convolution of diquark, valence quark and 
sea quark distributions u(x,n) in the incident particles and the fragmentation functions 
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G(z) of quarks and diquarks into secondary hadrons. The diquark and quark distribution 
functions depend on the number n of cut Pomerons in the considered diagram. 




Figure 1: Cylinder diagram (cylinder is shown by dash-dotted curves) corresponding to 
the one-Pomeron exchange contribution to elastic pp scattering (a) and its cut which 
determines the contribution to inelastic pp cross section (b). Quarks are shown by solid 
curves and string junctions [20] by dashed curves. 

In the case of a nucleon target the inclusive spectrum of a secondary hadron h has 
the form pQ: 

= £ , (1) 

where the functions determine the contribution of diagrams with n cut Pomerons 

and w n is the probability of this process [2TJ. Here we neglect the contributions of 
diffraction dissociation processes which are comparatively small in most of the processes 
considered below. This diffraction dissociation contribution is important mainly for sec- 
ondary production in large Xp region, which in not important in the present calculations. 
For pp collisions 

= n)f!(x-,n) + /*(*+, n)/* n) + 2(n - n)tf(x-, n) , (2) 



x± = - Am\j s + x 2 ± x] , (3) 

where f qq , f q and f s correspond to the contributions of diquarks, valence quarks and sea 
quarks, respectively. 

They are determined by the convolution of the diquark and quark distributions with 
the fragmentation functions, for example, 

f£{x+,n)= [ u q (x 1 ,n)G h q (x+/x 1 )dx 1 . (4) 
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The diquark and quark distributions as well as the fragmentation functions are de- 
termined from Regge intercepts [9]. 

In the case of nuclear targets we should consider the possibility of one or several 
Pomeron cuts in each of the v blobs of hadron-nucleon inelastic interactions as well 
as cuts between Pomerons. For example, for a pA collision one of the cut Pomerons 
links a diquark and a valence quark of the projectile proton with a valence quark and 
diquark of one target nucleon. Other Pomerons link the sea quark-antiquark pairs of 
the projective proton with diquarks and valence quarks of another target nucleons and 
with sea quark-antiquark pairs of the target. 

For example, one of the diagram for inelastic interaction with two target nucleons 
is shown in Fig. 2. In the blob of the pN\ inelastic interaction one Pomeron is cut, 
and in the blob of pN2 interaction two Pomerons are cut. It is essential to take into 
account every possible Pomeron configuration and permutation in all digrams. The 
process shown in Fig. 2 satisfies the condition |22j— [25j that the absorptive parts of 
hadron-nucleus amplitude are determined by the combinations of the absorptive parts 
of hadron-nucleon interactions. 
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Figure 2: One of the diagrams for inelastic interaction of an incident proton with two 
target nucleons N\ and N2 in a collision. 

In the case of inelastic interactions with v target nucleons n be the total number 
of cut Pomerons in hA collisions (n > v) and let rij be the number of cut Pomerons 
connecting with the i-th target nucleon (1 < < n — v + 1). We define the relative 
weight of the contribution with n« cut Pomerons in every hN blob as w\ . For the 
inclusive spectrum of the secondary hadron h produced in a pA collision we obtain |5j 
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j A f oo n—u+1 n—v+1 u 

= E v pa E E ■•■ E n< >< (5) 

(T pA "J-F v= x y n =v m=l n„=l Z=l 

x [ft q (x+, n)fi(x-, m) + n)/£(x_, nj) + 



2n-2 



+ E fs( x +i n )ft, q A x -> n ™)} } 



m=l 



where V^2 is the probability of "pure inelastic" (nondiffractive) interactions with v target 
nucleons, and we should account for all possible Pomeron permutation and the difference 
in quark content of the protons and neutrons in the target. 

In particular, the contribution of the diagram in Fig. 2 to the inclusive spectrum is 



%e da 
^fdx~F 



= 2V$< Nl wf 2 {/£(*+, 3)/ s fc (x_, 1) + (6) 

+ fi(x + , 3)/£(x_, 1) + f?(x + , 3)[/* (x_, 2) + f^(x., 2) + 
+ 2/ s ft (x_,2)]} . 

The diquark and quark distributions as well as the fragmentation functions are here 
the same as in the case of a nucleon target. 

3 Inclusive spectra in heavy ion collisions at inter- 
mediate energies 

At comparatively low energies we do not have a simple model for the calculation of 
the yields of secondaries. However, there exist the data [151 US], which show that the 
yields of secondaries produced by one interacting nucleon increase with the increase of 
the number of interacting nucleons. An example is presented in Fig. 3 taken from [16]. 

It is clear that the increase of the number of projectile participant, N part , means the 
decrease of the impact parameter and every participant interacts, on the average, with 
larger number of target nucleons, una [26J. One can see from Fig. 3 that every projectile 
nucleon produces more both K + and K~ at small impact parameters in comparison with 
peripheral interactions. These ratios change about 2-2.5 time from large to small impact 
parameters^. It means that the multiplicity of the produced kaons is a real function of 
vabi that is the number of elementary nucleon-nucleon collisions in the case of inelastic 
interaction of nucleus A with nucleus B and 

Kv AB ) - . (7) 

® AB 



1 Due to the energy correction factors [22] these ratios should not be equal exactly to vna 
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Fig. 3. The total yields of kaons per projectile participant versus the number of 
projectile participants N part in Au — Au reactions at 11. 6A GeV/c. 



This fact confirms the applicability of the Multiple Scattering Theory at the energies of 
GSI, as a minimum. 

The similar situation probably takes place at CERN SPS energies, see Fig. 10 in |17j . 
The saturation of the yields of secondaries produced by one interacting nucleon possibly 
can be seen only at N pp > 300, i.e. in very central events. However, it is necessary to 
note that the data [TTj are in some disagreement with the data [18j. 

At CERN SPS energies 158A GeV/c, following to the data p2], we can use the QGSM 
and calculate the spectra of secondaries. In the case of heavy ion collisions the energy 
conservation effects violate the asymptotical ratio 

1 da(AB -> hX) 1 da(NN -> hX) 

~p77d ~J7. - \ Vab I77m 777, 

a AB dy a NN dy 

significantly stronger than in the case of hadron-nucleus collisions, where the asymptot- 
ical value is 

1 da(NA -> hX) 1 da(NN-^hX) 

~^77d -J7, \ u na)— -r (9) 

<r NA dy a NN dy 
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with 



4 mei 

M = ■ do) 

The possible way to account for the energy conservation effects was suggested in [27J. 
The idea is that we will use the rigid target approximation in two different ways. In 
the case of forward hemisphere (nucleus A fragmentation region) we account for that 
each nucleon of the nucleus A can interact with several nucleons of the nucleus B, but 
the nucleons of B interact not more than once. This is equivalent to the case, where A 
uncoupled nucleons, with the corresponding impact parameter distribution, interact with 
the nucleus B. In the backward hemisphere (nucleus B fragmentation region) we will 
use the same, but with change A and B. Now B uncoupled nucleons interact with the 
nucleus A. The two contributions in the central region are matched with good accuracy. 

The inclusive spectrum of secondaries produced in A — B collision can be written as 



+ 9{-y)R h B {-y)(N B )- 



da{NA -> hX) 



a 



prod 
NA 



dy 



where y is the rapidity of secondary h in cm. frame and the functions R\ B {y) account 
for the energy conservation effects. 

Here we connect the inclusive spectra of secondaries in the heavy ion collisions with 
the spectra in the nucleon-nucleus interactions. To calculate the last ones, as well as the 
functions R\ B (y) we use QGSM and the Multiple Scattering Theory for NA collisions 
0. 

The functions R\ B (y) can be taken in the form [27] 

RAy) - ' (12) 

where f h {y,v) is the contribution to the secondary h spectrum from a beam nucleon 
interaction with v target nucleons, and f h (y, 1) is the spectrum of secondary particle h 
in NN collision. 

The rapidity distributions of secondary tc~ , K + and K~ [29J, as well as 7r + , p and p 
[T7j measured in Pb — Pb central collisions at 158 GeV/c per nucleon are compared with 
our calculations using Eq. (11) in Fig. 4. 

In the case of pions the agreement is reasonable. The differences between the calcu- 
lated curves and the data are not larger than 10%. In the case of secondary kaons we see 
that there is a problem. The spectrum of K~ also is in reasonable agreement with the 
model, whereas for K + disagreement is about 30%. Even more serious is the K + / K~ 
ratio. Experimentally this ratio at y ~ is about 1.5, however in the model we can not 
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Fig. 4. The rapidity distributions of secondary n + , K + , K~ , p and p produced in 
the central Pb — Pb collisions at 158 GeV/c per nucleon. Triangles show the results for 
the positive secondaries and over turned triangles for the negative ones. The solid 
curves are the QGSM predictions for positive secondaries and dashed curves for the 

negative ones. 



obtain more than 1.2, the ratio about 1.5 being obtained for pp collisions at the discussed 
energy. Most probably it means that the difference between fragmentation functions of 
diquarks and/or quarks into K + and K~ is not large enough in the model. The used 
fragmentation functions were taken from [2J. In the cases of p and p the agreement with 
the data is good enough. The contribution from string junction diffusion [30]- [33] to the 
proton spectrum is small at this energy. 

4 Inclusive spectra at high energies 

The Multiple Scattering Theory allows one to obtain some simple and model independent 
formulae [MJ [35] coming only from the assumption that high energy heavy ion collision 
can be consider as the superposition of independent nucleon-nucleon collisions. So any 
serious disagreement of these predictions with the data can be considered as the signal 
for some collective interaction. 



dN/dy 




Pb-Pb -» K 
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Fig. 5. Diagrams for inclusive cross sections for A — B collisions in Glauber 
approximation (a) and with accounting for the interactions of pomerons (shown by 
wave curves) (b). An example of inelastic processes of pp interactions which determine 
one of the vertices of the multipomeron interactions (c). 



The predictions for rapidity spectra of secondaries are more model dependent. The 
RHIC experimental data give clear evidences for the effects of inclusive density satu- 
ration, which reduce the inclusive density about two times in the central (midrapidity) 
region in comparison, say, with predictions of [36J, [37] based on the superposition picture. 
The observed phenomena can be explained in the framework of the space-time picture 
of high energy interactions and on Gribov's Reggeon diagram technique [SB]- In high 
energy hadron-nucleus collision there exists inelastic screening corrections [SHI SO]- The 
same inelastic screening should exist in high energy heavy ion collision. This effect is very 
small for integrated cross sections (because many of them are determined by geometry), 
but it is very important [19] for the calculations of secondary multiplicities and inclusive 
densities. 

At not very high energies the heavy ion collisions can be considered with the help 
of standard Glauber approximation, and the inclusive spectrum of any secondary h 
produced in the central region is described by the diagram shown in Fig. 5a, i.e. by the 
contribution of a single nucleon-nucleon blob. This diagram immediately gives Eq. (11) 
(and Eq. (12) after accounting the corrections for finite energy) for heavy ion inclusive 
cross section. The contributions of all other diagrams cancel each other due to AGK 
cutting rules [8]. 

At high energies the new diagrams appear, which include the interactions of pomerons 
and correspond to the diffractive production of a large mass M jet in pp interactions. 
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At low energies the contribution of such diagrams to nucleus-nucleus interaction is sup- 
pressed by the longitudinal part of nuclear form factor G z (t m i n ) which is connected with 
the longitudinal part of momentum transfer, q z (t m i n = —q^). Thiz longitudinal form 
factor can be written for Gaussian distribution of nuclear density p(b, z) as 

G A {t miTl ) = J p A (b, z)e^ z dz w exp (R 2 A t min /3) . (13) 

When energy becomes high enough, the value of t min becomes very small. So the dis- 
cussed contribution can be significant. 

One example of a diagram with pomeron interaction for heavy ion interaction is 
shown in Fig. 5b. Contrary to the hadron-nucleus case, where the inelastic screening is 
connected with the diffractive dissociation of projectile particle [391110], the contribution 
of such diagram can be estimated from the processes of high mass jet production in 
midrapidity region and with two large rapidity gaps, see for example Fig. 5c. The 
contribution of the considered diagrams to inclusive spectrum is suppressed quadratically, 
by both longitudinal form factors, GU(£mm) and Gs(t m m)- So we can observe their 
influence at energies quadratically higher in comparison with the energies region, where 
the inelastic screening effects are observed in hadron-nucleus scattering. 

Following to the estimations of [19], the RHIC energies are of the needed order of 
magnitude. The inelastic screening can decrease [19] the inclusive density in the midra- 
pidity region about two times at RHIC energies and about three times at LHC energies 
in comparison with the calculation without inelastic screening. 

However all such estimation are model dependent. The numerical contribution of 
all multipomeron diagrams is rather unclear due to a lot of unknown vertices for the 
multipomeron interactions, and the number of multipomeron diagrams is very large. 
The number of parameters can be decreased in some model, for example, in [19] the 
Schwimmer model [JT] was used for the numerical estimations. 

Another (again model dependent) possibility to estimate the contribution of the di- 
agrams with Pomeron interaction comes pl2| H3| HI] from percolation theory. In this 
approach we assume that if two or several pomerons are overlapping, they become a 
one pomeron. When all Pomerons are overlapping, the inclusive density is saturated, 
it reaches its maximal value at given impact parameter. This approach has one free 
parameter - the critical number of pomerons in one squared fermi. Technically it is 
more simple to bound the maximal number of pomerons, n max , which can be emitted by 
one participating nucleon for the given pair of colliding nuclei. All model calculations 
become rather simple because above the critical value every additional pomeron cannot 
contribute to the inclusive spectrum. 

The results of calculations of the inclusive densities of the produced charged secon- 
daries per one pair on interacting nucleons dn c h/dy / (0.5 N part ), with the help of Eq. (11) 
and QGSM are shown in Fig. 6. We present the CERN SPS experimental data for 
Pb — Pb collisions at ^Js = 17.3 GeV per nucleon [35]- [IB] f° r \y\ < 1 an d RHIC data for 
Au — Au at y/s = 130 GeV per nucleon [IB! EH] (these data are in agreement with the 
results of more recent measurements) . 
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Fig. 6. Relative inclusive densities of secondaries for Pb — Pb collisions at y/s = 17.3 
GeV per nucleon (black points, multiplied by 1/2) and for Au — Au at y/s = 130 GeV 

per nucleon (open points). Dashed curves present the QGSM results without 
percolation effects for y/s = 17.3 GeV (curve 1, multiplied by 1/2, as the points) and 

for y/s = 130 GeV (curve 2). The solid curve shows the results of calculations at 
y/s = 130 GeV with percolation effects ((n max ) = 1.67). Dash-dotted curves 1 and 2 
show the predictions for Pb — Pb collisions at 5.5 TeV per nucleon with (n max ) = 1.67 

and (n max ) = 1, respectively. 



One can see that CERN SPS data at y/s = 17.3 GeV are described reasonably (dashed 
curve 1 in Fig. 6) without any additional screening corrections, by the same way as the 
results shown in Fig. 4 for different secondary hadrons. It seems to be natural because 
the suppression effects of t min , Eq. (13), discussed above should be very important at 
comparatively low energy of CERN SpS. 

However, the same calculation at RHIC energy y/s = 130 GeV per nucleon gives 
the relative inclusive density two times larger (dashed curve 2 in Fig. 6) than the data. 
The agreement with the data can be obtained only with suppression of multipomeron 
contributions, namely by using (n max ) = 1.67 for every interacting nucleons. In this case 
the averaged number of cutted Pomerons emitted by every interacting nucleon is about 
1.2 that is significantly smaller than the number of cutted Pomerons in pp collisions at 
the same energy (the last one is about 1.8). 

Here we assume that the value (n max ) does not depends on N part in the interval A/4 < 
N par t < A, i.e. between minimum bias and very central collisions [26]. This assumption, 
as one can see from Fig. 6, is in agreement with the presented RHIC experimental data 
with the accuracy 10-20 %. It is confirmed (within the same accuracy) by the recent 
data [50] that the ratio of charged particle multiplicities in central (where N part « A) 
Au — Au and Cu — Cu collisions in midrapidity region are both at y/s = 200 GeV and 
62.4 GeV only 15 % large that the ratio of Au and Cu atomic weights. However, it is 
necessary to note that the calculated inclusive densities in heavy ion collisions per one 
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pair of participants are about 1.6 times larger than in pp collisions. 

The same value of (n max ) = 1.67 allows us [51] to describe the PHOBOS point at 
y/s = 56 GeV [H] . It means, that the saturation effects at y/s = 56 GeV are of the same 
order of magnitude as at y/s = 130 GeV and 200 GeV. 

The result of PHOBOS Coll. [52j for central Au - Au collisions at -y/s = 200 GeV 
per nucleon gives dn c h/drj = 650 ± 35 for \r)\ < 1 that is in agreement with our result 
dn c h/dr] m 600. The important point is that the experimental data [S3] of NA50 Coll. for 
central Pb — Pb collisions show the increase of the multiplicities of secondaries, dn c h/dr], 
in the energy interval y/s = 8.8 GeV - 17.3 GeV about two times, from 207 ± 1 ± 16 
to 428 ± 1 ± 34. In the much larger (in logarithmical scale) interval y/s = 17.3 GeV - 
200 GeV these multiplicities increase only about 1.5 times. Probably it means that the 
inelastic screening (percolation) effects start to work at energies about y/s ~ 20 — 30 
GeV and they become very significant at y/s ~ 50 GeV. 

More detailed data are presented in [51] and in [55], where the multiplicities of iden- 
tified hadrons were measured in Au — Au at y/s = 130 GeV and 200 GeV, respectively. 
The values of dn/dy/ (0.5 N part ) for different secondaries produced in 5 % of the most 
central Au — Au interaction are compared with our QGSM calculations in the Table. 
Percolation effects were accounted for, as it was explained above. The agreement of the 
QGSM calculations with the data is on the level of 20 % that is usual for QGSM (let us 
note that we did not input any new parameter). In particular, one can see qualitative 
agreement of the calculated results and the data in the energy dependences of midra- 
pidity multiplicities for secondary pions and kaons, but rather strange disagreement for 
secondary protons and antiprotons. 

Table 

The values of dn/dy / (0.5 N part ) for different secondaries produced in 5 % of the most 
central Au — Au collisions at y/s = 130 GeV and 200 GeV. 



Hadron 


130 GeV 


200 GeV 






54 




QGSM 




55 




QGSM 


7T + 


1.59 ±0.05 


1.82 


1.63 ±0.13 


2.00 




1.55 ±0.05 


1.88 


1.61 ±0.13 


2.05 


K + 


0.27 ±0.02 


0.17 


0.28 


±0.04 


0.19 


K~ 


0.23 ±0.02 


0.17 


0.26 ±0.04 


0.18 


P 


0.16 ±0.01 


0.13 


0.10 ±0.01 


0.14 


P 


0.11 ±0.01 


0.08 


0.08 


±0.01 


0.10 



The results of Dual String Model [561 EZ] calculations with string fusion [51] are in 
agreement with the results of the percolation model. The results [58] of String Fusion 
Model [59] are also in agreement with presented calculations. 
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Fig. 7. Total number of charged particles detected within range —5.4 < r/ < 5.4 in 
Au — Au collisions at ^/s =130 GeV per nucleon as a function of N part and its 
description by the QGSM with percolation effects and {n max ) = 1.67 solid curve. 



In the case of LHC energy a/s = 5.5 TeV per nucleon for Pb — Pb collisions with 
(n ma x) = 1-67 percolation effect decreases the relative inclusive density about 3 times. 
This result is shown by dash-dotted curve 1 in Fig. 6 and it is again similar to [19] 
prediction. Let us note that the energy difference in calculations with and without 
percolation effects is connected mainly with the increase of the total inelastic NN cross 
section. 

However, the percolation effect at LHC energy can be even larger if the increase of 
the transverse range of a pomeron with incident energy will be accounted for [51]. The 
predictions for the relative inclusive density at LHC energy with (n max ) = 1 (maximal 
percolation) are shown in Fig. 6 by dash-dotted curve 2. Now the percolation effect 
decreases the relative inclusive density about 5 times. 

The data [60] on the total number of charged particles detected within range —5.4 < 
7] < 5.4, i.e. in the region which exclude only high-x^ fragmentation and diffraction 
regions are presented in Fig. 7. 

These data are in good agreement with the QGSM calculations accounting for the 
percolation effects. The calculated dependence are very close to a straight line that is 
the direct consequence of Eq. (11). 

5 Conclusion 

We can see that the QGSM together with Multiple Scattering Theory can describe on 
reasonable level the inclusive spectra of different secondaries produced in heavy ion 
collisions at not very high energies. Some disagreement in K + and K~ yields in Fig. 4 
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can be corrected by reasonable change of the diquark and quark fragmentation functions 
into different kaons. 

The data of RHIC and their comparison with CERN SPS data show numerically 
large effects coming from the Pomeron (secondary particle) density saturation. It is the 
first experimental evidence of so numerically large effects in high energy physics. The 
inelastic screening effects should be accounted for in calculations of inclusive spectra of 
different secondaries produced in heavy ion collisions at RHIC and LHC energies. The 
processes of baryon number transfer via string junction diffusion [30j-[33j also should be 
accounted at these energies. 

From the comparison of the RHIC data (where effects are large) and CERN SPS 
data (where effects are small) we can conclude that the inelastic screening (saturation 
effects) become very important for heavy ion collisions in the energy interval ^/snn = 
20 - 30 GeV and they, probably are saturated at ^/snn = 60 GeV. At higher energies 
the growth of total inelastic NN cross section results in the increase of the number 
of Pomerons which can be screened, so the difference in calculations with and without 
inelastic shadowing becomes larger. 

We are grateful to N. Armesto, A. Capella, C. Pajares and R. Ugoccioni for useful 
discussions. This paper was supported, in part, by grants RSGSS-1124.2003.2 and PDD 
(CP) PST.CLG980287. 
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